Demonstration of superluminal effects in an absorptionless, non-reflective system 
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We present an experimental and theoretical study of a simple, passive system consisting of a 
birefringent, two-dimensional photonic crystal and a polarizer in series, and show that superlu- 
. . . minal dispersive effects can arise even though no incident radiation is absorbed or reflected. We 

CO ' demonstrate that a vector formulation of the Kramers-Kronig dispersion relations facilitates an 

' understanding of these counter-intuitive effects. 
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(N . 

^jpj' Superluminal group velocities have been observed in a number of different physical systems. These include passive 
^ \ absorptive , passive reflective [2, ll| , and active transparent Q media. There have also been numerous theoretical 
.^fH . and experimental proposals to observe superluminality in the tunneling of electromagnetic wavepackets Q. 
^ ■ Here we report the first experimental observation of superluminal effects in a passive system with neither absorption 
,; nor reflection. The effects arise because of a transfer of energy or interference between two modes of the electromagnetic 

■ field, in this case two different polarizations of light. We are able to interpret these results using a new vector 
^ , formulation of the Kramers-Kronig relations. 

' It is widely believed that the Kramers-Kronig (K-K) relations |^ require that passive systems be either absorptive 
ly-j \ or reflective in order to exhibit superluminal effects. In this paper, we demonstrate that this is not the case; while 
CO ■ absorptive and reflective systems are required to have spectral regions of anomalous dispersion, the converse is not 
T-H ' necessarily true. In fact, the exchange of energy between modes is a sufficient condition for superluminal propagation 
in any system. We show that these effects are consistent with causality. 
' Our experimental system consists of a slab of highly birefringent two-dimensional (2D) photonic crystal and a linear 
polarizer, placed in series. The photonic crystal has fundamental and second-order photonic band gaps in the regions 
of 10 and 20 GHz, and displays strong birefringence with very high transmission in the frequency range between the 
two gaps (7(1 . The crystal itself is an 18-layer hexagonal array of hollow acrylic rods (outer diameter 1/2") with an 
air-filling fraction (AFF) of 0.60. The crystal was constructed using a method which we have previously described j^. 

We studied the transmission and dispersive properties of this system between the two band gaps, using an HP 
8720 A vector network analyzer (VNA). Microwaves were coupled to and from free space with polarization-sensitive 
^2 ' horn antennae. The photonic crystal was placed in the far-field of the transmitter horn to ensure that planar wavefronts 
of a well-defined polarization were incident on it. The receiver horn was positioned immediately behind the photonic 
IL* crystal on a direct line of sight with the transmitter horn. In addition, the crystal and receiver horn were placed inside 
. !^ ] a microwave-shielded box with an open square aperture, whose size was chosen to minimize diffraction effects while 

■ eliminating signal leakage around the crystal. This method has proven very effective for transmission measurements 
?H ' at centimeter wavelengths 

In order to control the polarizations of the incident and detected fields relative to the fast axis of the crystal, we 
mounted the transmitter and receiver horns on precision rotation stages. The angle of the incident polarization 9 
was held fixed at 45° relative to the fast axis of the crystal, while the angle of the receiver horn j3 was allowed to 
vary. We define our coordinates such that the slab is oriented with its fast axis parallel to the vertical direction, and 
label the incident polarizations TM (transverse magnetic) and TE (transverse electric) for polarizations parallel and 
perpendicular to the fast axis, respectively. Our experimental setup is displayed in Fig. 

Our transmission amplitude and phase measurements with the receiver horn set at 40° and 50° are shown in Figs. 
|2Ia) and|2Ib), respectively. Since this spectral region is far from the band gaps, the transmission dip is not due to 
any band gap effect. It is caused by the fact that the photonic crystal rotates the polarization of the light by adding 
a different frequency-dependent phase to each polarization component 0- For (3 < 45° there is clear anomalous 
dispersion in the vicinity of 16.5-17 GHz (the half-waveplate frequency for the photonic crystal) while at /3 > 45° the 
dispersion is normal. 

A remarkable feature of our data is that while the transmission is identical for both /3 — 40° and /3 — 50°, aside 
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FIG. 1: Experimental setup. 
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FIG. 2: Calculated (dotted lines) and measured (solid lines) transmission (left axis) and phase (right axis) for the detected 
polarization inclined at (a) 40° and (b) 50°. 



from experimental errors, the phase properties are quite different. Motivated by these pecuhar results, we would first 
like to derive a simple physical model capable of explaining these phenomena and illustrating their connection with 
causality. 

A linear system which is invariant under time translation and is described by a scalar response (Green's) function 
g{t) produces a time-dependent response b{t) to an input a{t) given by the convolution b{t) = f^°° g{t — T)a{T)dT. 
This expression is the starting point for the usual derivation of the Kramers-Kronig relations ^ . Our system, however, 
accepts a time-dependent vector input and produces a vector output. We must therefore replace the scalar convolution 
with the expression 
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/oo poo 
gij{t - T)aj{T)dT ^ I gij{t')aj{t-t')dt' (1) 
oo J — oo 

where hi{t) and ai{t) are the i*'' components of the time-dependent output and input vectors b(i) and a(t), gij{t) is 
the {i,jY^ component of the Green's function matrix which describes the system. We have employed the standard 
summation convention over repeated indices. In the frequency domain, Eq. Stakes the simple form 



B,M = G,,Hi,H (2) 

where Bi{uj), Aj{uj), and Gij{uj) are the Fourier transforms of bi(t), ai{t), and gij{t), respectively (assuming these 
transforms exist). 

Causality tells us that the output must vanish for times before the input has propagated through the system. Thus, 
we require j'^^ gijiT)aj{t — T)dT — where T must be greater than or equal to the relativistic propagation time. 
Since this relation must be satisfied for any choice of input, the components of the Green's function matrix vanish 
individually for all time prior to T . This directly implies that the components (u;)e~*""^ have analytic continuations 
for complex frequencies u> where Imiu)) > (i.e., the upper half-plane in complex frequency space). Therefore, it is 
possible to proceed in the usual way to show that the real and imaginary parts of each component Gij{uj) satisfy the 
Kramers- Kronig dispersion relations assuming the components Gij{uj) are square integrable: 

ReC.,M^ipr%*2p« (3) 

TT Jo Vl^ ~UJ^ 



where P denotes Cauchy's principal value 0. 

The signals in our experiment can be expressed as A(lj) = A{Lj)eA and B(lj) = B[uj)eB, where and are unit 
vectors in real space, and ^(i^) and B(lu) are the complex amplitudes of the incident and detected electric fields. In 
this special case where and are frequency-independent, Eq. |21 implies that a scalar transfer function describes 
the relationship between the complex functions B{lu) and A{(jj) 



Hico)^J-^=G,,ei,e], (5) 
A{uj) 

where e\ and are the i*'' components of the unit vectors ba and e^, respectively. Since the elements Gij{oj) 
satisfy the Kramers-Kronig relations, it follows that the scalar function iJ(w) must also satisfy them. However, the 
components of (w) interfere with each other in H{lj). 

It is this interference which leads to the superluminal effects in our experiment. The Green's function matrix for 
our (birefringent) system is 



where (pTsi^) — nTE{^)'-^d/c and (/)ta/(w) = nTM{^)^d/ c are the frequency-dependent phases imparted to TE 
and TM polarizations, riTEi^) and nxAiiuj) represent the frequency-dependent indices of refraction for the two 
polarizations, d is the slab thickness, and c is the vacuum speed of light. Substituting these matrix elements into Eq. 
|S1 we find 

H{u;) = sin(/3) sin(6')e"^^^('^) + cos(/3) cos(6l)e''^^^-^("). (7) 

For incident linear polarization at an angle of 6 — t:/ A from the vertical, the magnitude and phase of the transfer 
function are, respectively. 
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{1 + sin(2/3) cos [A(/)(w)]} 



1/2 



(8) 



arg[H(w)] = arg 



where A0(cj) = (j)TE{^) — 4>tm{^)- 



arctan 



sin[A(/)(tj)] 



i[^<p{Lo)] + cot(/3) 



(9) 



The transmission 



H{uj) has a minimum for A0(wm) = (2m + l)7r where is a half-waveplate frequency of the 

= since the polarization which emerges from the slab 



photonic crystal and m is an integer. For /3 — 7r/4, H(ll!^ 

has zero projection along the unit vector e^. At this point, the phase arg 
vanishes. For f3 — n/Azt e where < e <C 1, we find 



is undefined since the transmission 
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garg[g(u;)] 
dui 



± 



1 aA0(cj) 

2i duj 



(11) 



where Tg = d arg H{lo) / duj is the group delay of the transmitted wave. As e — > (i.e., as we approach the singularity 

in the transfer function), the transmission goes to zero at uj,n and the group delay is unbounded. For example, if 
(3 = 7r/4 — e and 9A(/)(ci;)/9a'|^ > 0, it is clear that the group delay can become superluminal, zero, or even negative 
depending on the value of e. Thus, if one were to measure the time of flight of an analytic pulse based on the arrival 
of its peak, superluminal results can be obtained. However, there is no violation of causality here: the group velocity 
of the pulse has nothing to do with the signal velocity (i.e., the velocity of "information"), the quantity restricted by 
relativity. As we will demonstrate empirically, these superluminal group velocities are in fact required by causality 
through the Kramers-Kronig relations. 

In Figs. Ufa) and|2fb), we have plotted the results of this simple model described by Eqs. |Sland|nifor comparison with 
the experimental results. The frequency-dependent phases 4>te{'-^) and 4>tm(!^) were calculated using an indepedent 
measurement of the indices of refraction of the photonic crystal, and include phase delays associated with free- 
propagation in the air-spaces between the horns and the crystal. It is clear that the model correctly predicts the 
behavior of this system. Some discrepancies between the model and the experiment are observed at the edges of the 
plotted regions due to the effects of the band structure of the photonic crystal, which were not included in the model. 

In order to transform between the real and the imaginary parts of H{uj), slightly modified versions of the Kramers- 



Kronig integrals must be used since the function 



is not square integrable over all real frequencies. In practice. 



however, one only has actual data for H{lo) over a finite range of positive frequencies between values we label as ui 
and ^2- In this case, it is possible to apply approximate transforms by truncating the integrals of Eqs. |31and^at oji 
and UJ2- These transformations are approximately correct for H(uj) given in Eq. \7\if uJi ^ lo <^ uj2. 



On the other hand, it seems counter-intuitive, from a Kramers-Kronig point of view, that the function 



IS 



symmetric around (3 — 7r/4, while argH(uj) lacks symmetry about this point. In particular, it is usually possible to 
obtain the phase of a causal physical response function given its amplitude using the equation [To| 



~ , , 2cj 
argiJ(a;) = (Iquj 



In H{n) 



(12) 



where do is a fixed, undetermined constant which characterizes the system in spectral regions of constant transmission. 
Clearly, this transformation cannot produce two different phase functions given a single magnitude. However, it is 
known that this expression must be modified if the function H{Cj) has zeros in the upper half-plane Thus, the 
phase of the response function is not uniquely determined from a measurement of the amplitude alone; if the response 
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FIG. 3: Phase data (dotted line) and the result of the amplitude-phase Kramers-Kronig relation applied to the transmission 
data (solid line). 



function has zeros in the upper-half plane, these zeros increase the group delay in a nontrivial way over what is 



expected simply from 



For the transfer function given in Eq. [7| we find that the zeros satisfy 

A0(w„) = -iln|cot(/5)| + 27r(n+ i) (13) 



where n is an integer and < (3 < 7r/2. Given our form of A(j){u!), it follows from Eq. ^|that all the relevant zeros 
lie in the lower-half plane for < (3 < n/A, while all are in the upper-half plane for tt/A < (3 < 7r/2. Under these 
conditions, we see that it should be possible to apply a relatively simple Kramers-Kronig transformation to obtain 
the phase of the response function given its magnitude if < /? < 7r/4, but not if 7r/4 < fi < 7r/2. 

In Fig. 131 we display the Kramers-Kronig transformation of Eq. 1121 applied to the amplitude data shown in Fig. 
|2Ia) with do = 1.34d/c (do was determined by comparing the result of the integral transformation with the actual 
phase data far from the half-waveplate frequency). The actual phase data are also shown in this plot for comparison. 
Clearly, the transformation works quite well at /3 ~ 7r/4 — 0.087, confirming that the transfer function does not have 
any complicating zeros in the upper half-plane at this angle. However, this simple transformation cannot produce the 
radically different phase data obtained for (3 ~ tt/A -\- 0.087 (shown in Fig. [2Ib)). 

The interference between terms in H{lj) arises because of the two (polarization) modes available to light in this 
system. The birefringent photonic crystal allows a coupling between these modes, and energy can flow from one to 
the other. This feature is common to all systems that display superluminal effects. In absorptive media energy is 
scattered from the incident field into other directions by fluorescence, while in transparent active media two frequencies 
are present and the medium provides a coupling between them. In passive, reflective systems energy is transferred 
between the incoming and reflected waves. Since the vector formulation of the K-K relations presented above can 
apply to any pair of modes, any medium that displays this coupling between an incident mode and another mode 
should also display superluminal effects. This result extends the conclusion of Bolda et al., that the (scalar) K-K 
relations imply superluminal group velocities must exist in systems that are absorptive over some range of frequencies 

El- . 

In conclusion, we have shown that superluminal and even negative group velocities can exist as a result of inter- 
ference, instead of absorption or reflection. Moreover, it is evident from this discussion that causal superluminal 
propagation can potentially be observed in any system in which input energy can "escape" into other modes. Our 
results also demonstrate an example of an unusual application of the amplitude-phase Kramers-Kronig relations; for 
the system we have described, an infinitesimal adjustment of parameters radically affects the validity of a simple 
transformation between the amplitude and phase. 
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